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We theoretically show that the dipole force of an ensemble of quantum emitters embedded in a 
dielectric nanosphere can be exploited to achieve near-held optical levitation. The key ingredient is 
that the polarizability from the ensemble of embedded quantum emitters can be larger than the bulk 
polarizability of the sphere, thereby enabling the use of repulsive optical potentials and consequently 
the levitation using optical near-helds. In levitated cavity quantum optomechanics, this could be 
used to boost the single-photon coupling by combining larger polarizability to mass ratio, larger held 
gradients, and smaller cavity volumes while remaining in the resolved sideband regime and at room 
temperature. A case study is done with a nanodiamond containing a high-density of silicon-vacancy 
color centers that is optically levitated in the evanescent held of a tappered nano-hber and coupled 
to a high-hnesse microsphere cavity. 


I. INTRODUCTION 

Both a two level atom and a lump of dielectric ma¬ 
terial are polarizable objects whose motion can be ma¬ 
nipulated with optical light [1], a feature that has been 
exploited in the fields of ultracold atoms [2] and opti¬ 
cal tweezers [3]. A recent experiment [4] has shown that 
the force from quantum emitters embedded in a dielec¬ 
tric nanosphere can be observed-even in liquid. In this 
article we go a step further and argue that it is exper¬ 
imentally feasible to have a scenario where the optical 
dipole force due to the embedded quantum emitters is 
stronger that the bulk dipole force. This would allow to 
use repulsive optical potentials for dielectric nanospheres, 
by using the internal structure of the quantum emitters, 
to trap them in evanescent helds analogously to what 
can be done with ultracold atoms [5]. We propose to 
use this for dispersive cavity quantum optomechanics [6] 
with optically levitated nanospheres [7-13]. In addition 
we show that it might be possible to reach the so-called 
strong single-photon cooperativity regime [6, 14-19], in 
the resolved sideband regime [6], by coupling the levi¬ 
tated nanospheres to a high-finesse optical micro cavity, 
at room temperature. 

In dispersive cavity quantum optomechanics [6], the 
single photon coupling go is proportional to 
where a is the real part of the polarizibility, b is the 
gradient of the cavity field mode, 14 is the cavity vol¬ 
ume and M the mass of the mechanical oscillator. Our 
proposal aims at simultaneously combining salient fea¬ 
tures that have been demonstrated in independent exper¬ 
iments: i) placing a dielectric nanobject in the near-held 
of a microcavity to have larger b and smaller 14 [20], 
ii) use a set of N quantum emitters to have larger 
a/Mi/2 oc ^/]V [21, 22], and hi) optically levitate a 
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nanosphere in high vacuum to have a high mechani¬ 
cal quality factor at room temperature [13] with a suf- 
hciently large trap frequency to enable cavity cooling 
in the resolved sideband regime [11]. This conjunction 
of features could be used for observing non-Gaussian 
physics of mechanical nanooscillators [14-19], measur¬ 
ing short-distance forces [23] , and migrating to quantum 
optomechanics cutting-edge experiments and proposals 
done with ultracold atoms in near-helds [24-33]. We re¬ 
mark that the increase of the polarizability of diamond 
by doping it with a high-density of color centers might 
also be used to boost the optomechanical coupling in se¬ 
tups using clamped cantilevers, membranes, or photonic 
crystals made of diamond [34-36]. See [37] for similar 
ideas. 

The paper is organized as follows. In Sec. II we dis¬ 
cuss and compare the polarizability of quantum emitters 
to the bulk polarizability of the nanosphere. We show 
that for a particular type of color centers in diamond, 
the polarizability of the quantum emitters can overcome 
the bulk polarizability. Assuming this regime, we de¬ 
scribe in Sec. Ill the optical trap obtained using a bi- 
chromatic held supported by a nano-hber. In Sec. IV we 
show that by placing a micro-cavity in close proximity to 
the nanosphere it is possible to achieve both the resolved 
sideband and the strong single-photon quantum optome¬ 
chanical cooperativity regime. We conclude in Sec. V. 


II. QUANTUM EMITTER POLARIZABILITY 
A. General expression 

Let us compare the polarizability of a quantum emit¬ 
ter (called quantum polarizability hereafter) with the 
bulk polarizability of a dielectric nanosphere. In gen¬ 
eral terms, the time-averaged dipole force describing the 
interaction of a monochromatic held E(r, t) of frequency 
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uj with a particle of polarizability a is given by [ 1 ]: 

r(r)=ai?o(r)^^^, (1) 

where E^iv) = 2(|E(r, t)|^) (here (•) denotes time av¬ 
erage). In the case of a bulk dielectric nanosphere of 
refractive index n and radius R, the polarizability as is 
given in the point dipole approximation by 

— 1 



where eo is the vacuum permittivity and V = AttR?/Z. 
The validity of the point dipole approximation, usually 
assumed by the condition R ^-kcIw for a focused Gaus¬ 
sian beam, can be verified in the context of near-fields us¬ 
ing a multimodal decomposition for the near-field [38, 39] 
and the dielectric sphere. Using typical experimental pa¬ 
rameters one can show that the interaction can be very 
well approximated by the point dipole term. 

For a two-level state quantum emitter with tran¬ 
sition frequency wo, dipole moment d, Rabi fre¬ 
quency O = |(e|d • e|(?)| i?o(r)/Ii, spontaneous emis¬ 

sion in free space Fq = (Ploq/(ZcottHc^) (accounting 
for an orientational average for the dipole moment 
|(e|d • ejg)]^ = cP/Z), spontaneous emission inside the di¬ 
electric nanosphere F k, nFo [40], and transverse decay 
rate 7 = r/2 -|- 7c, where 7c accounts for the additional 
coherence decay (inhomogeneous broadening), the quan¬ 
tum polarizability a^ is given by [41] 


2Ad2r s 
3 fiU ^7 s -I- 1 


( 3 ) 


Here A = w — wq is the detuning and s = H^/[ 7 F (1 -|- 
A^/ 7 ^)] the saturation parameter. This description uses 
the rotating wave approximation, which is valid pro¬ 
vided jAj <C Wo, and the Born-Oppenheimer approx¬ 
imation, which is valid provided the motional dynam¬ 
ics are much slower than the electronic dynamics of 
the quantum emitter. Note that Uq can be maximized 
to aq = —dP/{ZhA) for the optimal detuning jAj = 
7[1 -I- H^/( 7 r)]^/^ that leads to a saturation parameter 
s = 0 ^(r 7 )“^(l -I- A^/ 7 ^)“^ < 1. The optimal ratio be¬ 
tween the two polarizabilities 77 = \aq/as\ can thus be 
written as 

^ 2 rP+2 T 1 

^ i?3 (4,^)3 rp -In'y 


Assuming that in the nanosphere one has a number of 
embedded identical quantum emitters given by A = PqV, 
where pq is the volume density of emitters, and that the 
variation of fields within the nanosphere is negligible such 
that each quantum emitter identically interacts with the 
field (point-particle approximation), the total quantum 
polarizability will be larger than the bulk polarizabil¬ 
ity when Nrj > 1 . Alternatively, an effective complex 


FIG. 1: Ratio Nr) (solid dark line) of the total quantum polar¬ 
izability over the dielectric polarizability for a nano-diamond 
with a high density of SiV. Larger intensities both increase the 
saturation of the quantum emitters and the internal temper¬ 
ature (solid light line) leading to a reduction of the quantum 
polarizability. The properties of the nanodiamond and the 
SiV are listed in Appendix A 1 and A 2 respectively. We re¬ 
mark that since the SiV centers are embedded in the ND, we 
accounted for the Lorentz local field correction [51], (n^-|-2)/3, 
for the field inside the ND. 


refractive index n can be defined from the complex po¬ 
larizability dq of the ensemble of emitters by using the 
Lorentz-Lorenz relation pqdq/{Zep) = {fp — l)/(n^ -I- 2) 
for which one would obtain a value close to 2i near the 
resonance of the quantum emitters. It is important to 
note that we have not taken into account cooperative ef¬ 
fects. These effects could have a significant impact on 
the dipole force and consequently polarisability ratio 77 . 


B. Color centers in diamond: silicon-vacancy 
centers 

Let us consider a levitated nanodiamond (ND) [43- 
45] for its remarkable optical properties and the large 
variety of color centers they have [46, 47], which act 
as quantum emitters. In particular, the silicon-vacancy 
(SiV) centers, consisting of a silicon atom and two ad¬ 
jacent vacancies, have the following favorable properties 
for achieving Nr) ^ 1: i) strong dipole moment [48], ii) 
high-densities [48] , and iii) very good indistinguishability 
of different SiV [49] . 

The polarizability of the quantum emitters varies lin¬ 
early with the intensity of the incident field for small 
intensities (H^ <C yF) leading to a constant ratio 77 . For 
larger intensity, this ratio decreases due to saturation ef¬ 
fects from the quantum emitters. The incident intensity 
also has an impact through the dependence of the emit¬ 
ters with the internal temperature of the ND [50] (see 
Appendix A 2). In the context of levitation, this effect is 
important as the ND can reach high temperatures even 
with low intensities, see Fig. 1. In order to account for 
this effect, the internal energy of the ND is obtained as a 
function of the absorption rate of the incident laser pho¬ 
tons, 77 , the absorption rate of the black body radiation 
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of the environment, 7 a, and the black body emission rate, 
7 e, is given by [52] 


mCmT^ = 7/ + / dw' [7a (w') - 7e (w', Ti)] Huj', (5) 


where Ti is the internal temperature of the ND, m its 
mass, and Cm its specific heat. For an incident laser of 
intensity I and frequency w, and an environmental tem¬ 
perature Te, these rates are given by 




li = 


la = 


le = 


-Im 


£(^) - 

e(w) -|- 2 

4 (uj'R/c)^ j-K 
exp (^Huj j 1 
4 ^ 

TT V C 


Im 


exp 


e(a.') - 1 


huj' \ 


£(w^) -|- 2 

eico') - 1 


e{uj') + 2_ 


( 6 ) 

(7) 

(8) 


The steady state temperature (Ti = 0) is determined 
as a function of the incident laser intensity (see Fig. 1) 
using the rate-balance equation for the internal energy 
(Eq. (5)). Accounting for the temperature dependence of 
the linewidth, lifetime and transition of the SiV centers, 
see App. A 2, the total quantum polarizability for a ND 
containing a large density of centers clearly overcomes 
the bulk polarizability for low intensities (see Fig. 1). 
The internal temperature will be largely modified by the 
non-radiative properties of the embedded quantum emit¬ 
ters which we have not accounted for here. In particular, 
the main contribution will be from the limited quantum 
efficiency of SiV centres and, to a lesser extend, from its 
emission in the phonon sidebands. In addition, it is im¬ 
portant to note that the experimental values used to es¬ 
timate the emissivity of the diamond (see Appendix A 1) 
are limited to the 2.5-6.5^m range leading to an underes¬ 
timated black body emission rate. The internal temper¬ 
ature of the nano-diamond could be further reduced by 
coating its surface with a small layer of glass [44] provid¬ 
ing much higher emissivity in the infra-red. By achieving 
much lower internal temperatures, the optical properties 
of the SiV centers will significantly improve. 


III. NEAR-FIELD OPTICAL TRAP 

In the following we assume that the quantum polariz¬ 
ability dominates Nr] ^ 1 and discuss in a general way 
how this can be used to levitate a dielectric nanosphere 
in evanescent fields. In analogy to bichromatic atomic 
optical traps [53-55], we consider an evanescent field in 
vacuum of the form: 

E(r,t) = £;i(r)eie-‘"'i* + F;2(r)e'2e-'‘"^*+‘^ + c.c. ( 9 ) 

The exponentially decaying electric field amplitudes are 
given by T?i( 2 )(r) = , where Ai( 2 ) is the 

field decay rate and Ef(x,y) the field amplitude at the 
dielectric surface. The polarization of each mode is given 


by ei] 2 ) and tp is the phase difference between the two 
modes at t = 0. This field can be obtained from a nano¬ 
fiber and is different from the cavity field which will pro¬ 
vide the back-action. We consider the symmetric driving 
uji = ajQ — A and cu 2 = wo + A with A > 0 and define the 
beating frequency <5 = a ;2 — = 2A. We are interested 

in the interaction of such field with a two-level quantum 
emitter when far-detuning or weak driving for the two 
modes is not assumed. This is an unusual scenario for 
most of the atomic trapping experiments with evanescent 
fields [5]; notwithstanding, this was studied before both 
experimentally and theoretically, see for instance [56-60] 
and reference therein. One can analytically calculate the 
time-averaged dipole force considering that the beating is 
much faster than the dynamics of the mechanical motion. 
This is done solving a recursive equation obtained with a 
Floquet analysis. As shown below, one encounters that 
the total dipole force is not simply the sum of the forces 
that each mode would exert in the absence of the other 
mode but a more involved expression containing mixing 
terms that reflect the intricate interplay between the red 
and blue driving fields. 

In particular, for a bichromatic field (Eq. (9)), the 
Hamiltonian of the system in the rotating frame defined 
by the unitary U(t) = exp [iwi le)(ej t] reads 


H = HA le)(el — Mgg • E(r, — H.c.. (10) 

Here deg = (ejdjg) and r is the position of the two-level 
system. Applying the rotating wave approximation (valid 
provided A ojq) this Hamiltonian can be written as 

H = -hA\e){e\ + ^^^^ [l + S(r)e-'^‘] d++H.c. (11) 
where we have defined 

2 

^l(2)(r) = ■ D(2)7fl(2)(r), (12) 

which is also assumed to be real, and 5(r) = 

D 2 (r)/Hi(r). Hereafter we omit the r dependence to ease 
the notation. Assuming that the internal electronic dy¬ 
namics of the two-level system are much faster than the 
motional dynamics (Born-Oppenheimer approximation), 
the force is given by 


E = -V 


hill 


(1 -t 


) 


((7+) -I- H.C., 


(13) 


where the expected values are calculated for the elec¬ 
tronic steady state. Following [57] the steady state so¬ 
lution can be expanded as a Fourier series in terms of 6 
(Floquet’s analysis) such that 


OO 

u{t) = (a-) = Y. (14) 

n— — oo 

CO 

v{t) = (d+) = ^ (15) 

n— — oo 

OO 

w{t) = {a^)= Y (16) 

n— — oo 
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The Fourier coefficients are obtained by solving the op¬ 
tical Bloch equations in the steady state, namely 

0 = -[y+ \5{n+\/2)]Un+i^ (Wn+'B.Wn+l) , 

VLi 

0 = - [7 + i(5(n - 1/2)] - i— (w„-I- Sw„_i) , 

0 = -{T+ inS)wn-iil.iivn + ^Vn+i) 

+iOi -I-— r(5„o. (17) 


Then, by further assuming that S is much larger that 
the motional frequency (secular approximation), one can 
time-average Eq. (13) for the steady state to obtain 


F = --— (uo-fwo)--- -{vi + u-i). (18) 


Note that the time-averaged population of the excited 
state Pe will be given by 

^0 “b 1 ifli r ] ^ / M ^'1n^ 

Pe= -^- = —[U 0 -VQ + ^{U-i-Vi)\ (19) 


In addition, Eq. (17) leads to the following recursive 
equation for Wn 



FIG. 2: Trapping potential for a 30 nm ND accounting for 
the dipolar forces acting on the quantum emitters and the 
ND along with gravity and the Casimir-Polder forces (see 
Appendix A3 for the experimental parameters). The trap¬ 
ping potential obtained using the dipole force for the bi- 
chromatic field using the second order in the Floquet analysis 
(solid line) shows an important deviation from considering 
the dipole force from the two evanescent fields independently 
(light solid line). For reference, the force without accounting 
for the Casimir-Polder forces is represented with the dashed 
line. 


^n'^n T ^n'^n+1 T ^n'^n — l — dm (^0) 


where 




b 


n 


Cn 

dn 


_ 40 ^( 7 -b m( 5 )(l-b^^) 

472-hi8n7^-b^2(i_4„2)’ 

2Enl 

27 -I- iS{2n + 1 ) ’ 

2Enl 

27 + iS{2n — 1 ) ’ 

EdnO- 


( 21 ) 

( 22 ) 

(23) 

(24) 


By obtaining the set of Wn on can readily calculate 
and Vn- This recursive equation can be exactly solved by 
fixing a cutoff Af > 0 such that w_\f+i = = 0. The 

value of Af is chosen such that the value of the physical 
quantities that are calculated converge. 

In the lowest order in the Floquet analysis, namely for 
Af = 0, one obtains: 


to prevent the nanosphere from adsorbing to the surface 
producing the near-fields. Consequently it is necessary 
to account for the Casimir-Polder force to fully capture 
the trapping potential. In order to simplify the calcula¬ 
tion, we have approximated the surface as a semi-infinite 
half space of dielectric material with dielectric constant 
e (w). The Casimir-Polder potential for an object of po¬ 
larizability a{uj) at a distance z from the surface takes 
the form [61, 62]; 

fc pOO poo 

^ k — g{x,k) / 2k^c^\ e(ix)k — g(x, k) 

k + g(x,k) \ x^ J e(ix)k + g(x, k) 

(26) 

where 

g{x,k) = //j [e(ia^) - 1] + (27) 


27 1 -b Si -b S2 


(25) 


We defined Si( 2 ) = Ml + (^/t)] ^ § 1 ( 2 ) = 

Vlogfli( 2 )(r), and fli( 2 )(r) = 2 jd^g • di( 2 )Ei( 2 )(r)|/A 
The population of the excited state of the quantum emit¬ 
ter is given by pe = [1 — (1 -b Si -b S2)“^]/2. The to¬ 
tal optical force exerted on the nanosphere can then 
be estimated by using Ft ~ NFq + Fg, with Fs(r) = 
Q;sEo(r)V£’o(r)/2 the optical dipole force due to the bulk 
polarizability. 

In order to achieve near-held levitation, the blue- 
detuned Held (A > 0) provides repulsive optical forces 


For the polarizability of the quantum emitters with a 
transition frequency ujq we use cnqipj) = 2d^uiQ[ih(ijjQ — 
0 ;^)]“^. For a dielectric sphere of refractive index n 
and volume V we use as{uj) = 3eoV(n'^ — l)/(n^ -b 2). 
Considering silica for the semi-inhnite half spaces, see 
Appendix and [63], one can numerically calculate the 
Casimir-Polder force. The total Casimir-Polder poten¬ 
tial for a nano diamond containing N quantum emitters 
is obtained as 


[/cp(z) = NU^q^iz) + [//P(z). (28) 

Finally, one can then estimate the total force exerted on 
the nanosphere by using Ft Ri NFq + Fcp -b Fs -b Fg. 
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FIG. 3: Schematic illustration of the general scenario that 
is considered: a dielectric nanosphere with a set of two-level 
quantum emitters is trapped using a bichromatic evanescent 
field at a distance z from a surface. In the inset of the figure 
the level structure of an individual quantum emitter of transi¬ 
tion frequency ljo and excited state linewidth 7 is illustrated. 
The transition is driven by a bichromatic field with symmetric 
red and blue detuning =pA > 7 , with Rabi frequencies fli and 
Q 2 , respectively. The cavity field mode with vacuum Rabi 
frequency flc is detuned by Ac > 7 . 


Here Fcp = —VC/cp is the Casimir-Polder force and Fg is 
the gravitational force along the z-axis. 

The most efficient near-held trap is obtained using the 
optimal detuning, |A| = 7 [ 1 -|-H^/( 7 r)]^/^, for the quan¬ 
tum emitters to maximize the polarizability ratio rj. The 
trapping potential for the particular case of a 30 nm ND 
with embedded SiV centers is shown in Fig. 2. The possi¬ 
bility to levitate the nanosphere using near-helds, which 
provides much stronger optical forces, allows the use of 
low intensities and thereby the reduction of heating ob¬ 
served when using a focused Gaussian beam [64]. The 
dephasing of the SiV, varying as the cube of the inter¬ 
nal temperature [50] , remains sufficiently low to obtain a 
total quantum polarizability much larger than the bulk 
polarizability. On a practical level, the trap depth is par¬ 
ticularly adapted to the use of mobile optical fiber traps 
as they provide a cooling of the center-of-mass motion 
down to 30 K with a relatively simple apparatus [65] . The 
experimental parameters are described in Appendix A3. 


IV. CAVITY OPTOMECHANICS 


A. Optomechanical coupling 


(d). The total Hamiltonian is then given by: 


H t'+t' I\/I ..fc _\ 

— = Wto' 0 -\- Wo |e)(e| -I- Wc® a H- - — [aa -I- o' cr j . 


(29) 

Here we used the rotating wave approximation, which is 
valid provided |Ac|,Hc -C wq, where Wc — wq = Ac- The 
vacuum Rabi frequency is given by flc{z) = 
with Oc = 2d[fiWc/{5Vc€o)Y^'^/fi- and Ac the decay rate 
of the cavity’s evanescent field. Here 14 is the cavity 
mode volume, z = Zzp{b''^+b) is the displacement operator 
along the z-axis from the equilibrium position zt, with Zzp 
the zero point motion, and Wt the mechanical frequency. 
The coefficient ^ accounts for the particular mode shape 
and polarization type of the cavity such that Ecav = 
^[hiWc/{ycto)Y^'^ is the vacuum field at the cavity/vacuum 
interface [ 66 ]. By further assuming |Ac| > 7 flc one 
can make a Schrieffer-Wolf transformation to obtain 


H 


=WtPb -I- 


“ 1 a 7 “ “ 


|e)(e| -I- Wcd'^d 


4A, 


(30) 


-cr a'a. 


Since ^1/{4,Ac) 7 , the atomic level shift due to the 

cavity field will be irrelevant and thus can be safely 
dropped. By Taylor expanding the vacuum Rabi fre¬ 
quency around the equilibrium position, one obtains 


— =wtb^b + wq |e)(ej -|- Wcd^d — go —d^d(6^ -I- b) (31) 

fi uj[) 

By considering that the single-photon coupling go is much 
smaller than the frequency of the internal dynamics one 
replaces the (Tz by the (cr^) = wo calculated in the steady 
state considering the bichromatic driving. Therefore we 
arrive at the final single-photon optomechanical coupling 
Hamiltonian 


— =wtb^b -I- Wo |e)(e| + Wcd^a — goa^aCb^ -I- b). (32) 

h 

where the single-photon optomechanical coupling due to 
N quantum emitters is given by: 

, q 2 

g^o = -N{2pe - l)f 

Note that the coupling due to the quantum emitters 
would vanish should the quantum emitters be totally sat¬ 
urated (pe = 1 / 2 ), leaving only the contribution due to 
the bulk polarizability from the nanosphere [ 20 ]. 


Let us now address the optomechanical coupling be¬ 
tween the levitated nanosphere and an optical cavity 
mode (see Fig. 3). This is achieved by placing the sphere 
in the evanescent held of an optical microcavity mode 
with frequency Wc and creation (annihilation) operator 


B. Resolved sideband regime 

In order to cool down the center-of-mass of the me¬ 
chanical mode into the ground state using the optome¬ 
chanical coupling to the cavity [67-69] one requires the 





















6 



FIG. 4: Optomechanical coupling for a 30 nm ND while 
in the resolved sideband regime. The best optomechanical 
coupling is obtained at low-intensity as it provides a small 
excited state population (minimising the cavity scattering 
losses, see Eq. 34) and low inhomogeneous broadening thanks 
to a moderate internal temperature (maximising the coupling, 
see Eq. 33). For the parameters see Appendix A 5. 

resolved sideband regime wt > k. In this context, the 
use of evanescent field provides much stronger gradient 
compared to focused Gaussian beams, allowing for larger 
mechanical frequencies. To verify that the resolved side¬ 
band regime can be reached, it is necessary to consider 
the impact that quantum emitters have on the optical 
quality factor through the scattering of cavity photons. 
For N emitters, this is given by 

= (34) 

Here, we assumed the saturation parameter from the 
cavity vacuum field is much smaller than 1, Scav = 
^^H^(r 7 )“^ [1 -b (A/y)]”^ <C 1. In order to reduce this 
effect, we used a far detuned cavity, taking advantage of 
the excellent optical property of fused silica in the near 
infrared (see Appendix A 4). The cavity loss rate is then 
obtained as k = 2Trc/{QXcav) + k,sc- 

By maintaining the ND away from the cavity surface, 
the impact of the embedded quantum emitter on the 
quality factor allows to maintain a strong optical trap 
in the resolved sideband regime with an optomechanical 
coupling comparable to the state of the art (see Fig. 4 
and Appendix A5). Upon reaching a low center-of-mass 
temperature, the coupling can be further increased by 
using lower intensities, i.e. smaller dephasing (by main¬ 
taining a low internal temperature). As illustrated in 
Fig. 4, the coupling can be varied by one order of magni¬ 
tude by changing the total intensity while remaining in 
the resolved sideband regime. Since the position of the 
nanodiamond is not changed, the optomechanical cou¬ 
pling decreases for large intensities due to the increase 
of both the internal temperature and the excited state 
population. 

Another advantage of such bi-chromatic evanescent 


O 



FIG. 5: Optomechanical coupling and ND position. By vary¬ 
ing the intensity ratio + U), the position of the 30 nm 

ND can be controlled in order to the reach the optimum op¬ 
tomechanical coupling. As the ND is brought close the cavity, 
the cavity scattering losses is greatly increased by the pres¬ 
ence of the quantum emitters. 

trap resides on the possibility to accurately control the 
position of the ND by simply varying the ratio of inten¬ 
sity between the two fields. Consequently, it is possible 
to further increase the optomechanical coupling by plac¬ 
ing the ND closer to the cavity. Fig. 5 illustrates the 
evolution of the optomechanical coupling as a function 
of the intensity ratio. For each point, the total inten¬ 
sity is the lowest one still maintaining the resolved side¬ 
band regime. As the intensity of the repulsive field is 
increased, the ND is pushed further away from the nano¬ 
fiber, i.e. closer to the cavity. This allows to increase 
the optomechanical coupling up until the cavity scatter¬ 
ing losses reach values comparable to the unloaded cavity 
losses. For shorter cavity-ND distances, the ratio of the 
optomechanical coupling to the total scattering losses is 
reduced. In addition, higher intensities are required to 
achieve the resolved sideband regime, which increases the 
population of the excited state Pe and hence reduces the 
coupling (Eq. (33)). 

C. Strong single-photon optomechanical 
cooperativity 

The challenging strong-single photon optomechanical 
coupling regime is achieved when > K,Tm [6], where 
K is the decay rate of the optical cavity mode and F^ 
is the decoherence of the center-of-mass mechanical mo¬ 
tion. The leading contribution to F^ in optically lev¬ 
itated nanospheres is the recoil heating due to scatter¬ 
ing of photons [8, 10, 70]. In our case there will be 
the additional and dominant contribution due to scat¬ 
tered photons from the quantum emitters. One can 
estimate that F^ = F^ -|- Y^i=i^rn^ where the con¬ 
tribution from the quantum emitters is [71, 72] F^ = 
A(2/5)((x>oZ2p/c)^i?|c with the scattering rate from a sin¬ 
gle quantum emitter being = Fpe, and the contri- 
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bution due to the bulk polarizability [8, 10, 70] = 

X]i=i(2/5)(wi22p/c)^i?|2 j with the scattering rate given 

by Rtc,i = \asf Ef{rt){u}i/c)^/{12TTeoh) (r* is the trap¬ 
ping position of the sphere). Although much weaker, 
the recoil due to the black body radiation from the 
nanosphere can be obtained from the black body emission 
rate Eq. 8 as F^^ = f dw'[(2/5)(a;'z2p/c)^7e(w')]. One 
can now define the single-photon optomechanical coop- 
erativity C = gg/(KFm). We show below that by taking 
advantage of the quantum polarizability from quantum 
emitters embedded in the nanosphere, it is possible to 
obtain the strong single-photon optomechanical cooper- 
ativity regime C > 1 while maintaining the resolved side¬ 
band regime. 

In order to reduce the impact the quantum emitters 
have on the cavity losses, large detunings can be used 
to trap the ND. This approach, inspired from the far- 
off-resonance trapping (FORT) for cold-atoms [73], relies 
on the use of a far detuned near-field trap to reduce the 
excited state population pe. In this case, one constrain 
remains: the force from the quantum emitters should 
remain larger than the bulk one in order to provide blue- 
detuned repulsive fields. Another consequence of this 
configuration is the important reduction of the overall 
forces. To compensate for this, we propose to use an ad¬ 
ditional blue-detuned field from the microsphere cavity. 
Such field plays the role of the red-detuned field from 
the nano-cavity in preventing the ND from falling on the 
cavity under the Casimir-Polder forces. By using a far- 
of-resonance blue-detuned field, one can achieve a much 
stronger field gradient, further improving the trapping ef¬ 
ficiency. Despite the added complexity of using trapping 
fields from both the nano-fiber and the cavity, this con¬ 
figuration allows for maintaining the resolved sideband 
regime while significantly reducing the population of the 
excited state by using a large detuning. In particular. 
Fig. 6 shows the trapping potential (see Appendix A 6 
for the experimental parameters). The 30 nm ND is 
then maintained in the near-field of the cavity allowing to 
reach both the resolved sideband regime and the single¬ 
photon strong cooperativity C > 1, see Appendix for the 
set of experimental parameters. 


V. CONCLUSIONS 

In this article we have shown that the polarizabil¬ 
ity given by a set of quantum emitters embedded in a 
nanosphere can overcome its bulk polarizability. This 
immediately allows to use repulsive optical forces and lev¬ 
itate a nanosphere using evanescent fields. We have then 
discussed how this can be used for levitated quantum 
optomechanics using an evanescent coupling of the cen¬ 
ter of mass of the nanosphere with a microcavity mode. 
Due to a large polarizability to mass ratio, large cavity 
field mode gradient, and a smaller cavity mode volume, 
one could achieve the strong single photon cooperativity 
in the resolved sideband regime. A case study has be 



FIG. 6: Trapping potential for a 30 nm ND using two blue- 
detuned fields. Two blue-detuned fields, one from the nano¬ 
fiber and from the cavity are used to trap the ND. The large 
detuning, |A| = 1 x 10^® Hz, allows to maintain a very small 
excited state population (see Appendix A 6 for the experimen¬ 
tal parameters). To provide confinement on the transversal 
dimensions, an additional weak red-detuned field from the 
nano-fiber can be used. 


done considering a nanodiamond with a high density of 
embedded SiV color centers at room temperature. Also, 
this work could be extended to a recently discovered color 
center in diamond, the germanium-vacancy center [76], as 
it provides significantly better optical properties at room 
temperatures. 

We remark that as suggested in the recent experiments 
reported in [4] , collective effects from the high-density of 
quantum emitters embedded in the nanodiamond seem to 
be relevant. In particular, these effects could significantly 
modify the dipole force acting on the quantum emitters 
but also the total radiative emission of the emitters. The 
latter is a key point since otherwise the imperfect quan¬ 
tum efficiency of single quantum emitters could heat up 
the nanodiamond to unfeasible temperatures. The the¬ 
oretical understanding and potential applications of col¬ 
lective effects in this scenario, namely a levitated sphere 
smaller than the optical wavelength with such a high- 
density of quantum emitters that dipole-dipole interac¬ 
tions have to be taken into account [74, 75] in addition 
to inhomogeneous broadening, is a very interesting and 
challenging further research direction that will be ad¬ 
dressed elsewhere. 
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Appendix A: Experimental parameters 


2. Quantum emitter: silicon-vacancy 


In the following we list all the experimental param¬ 
eters needed to make the analysis of the experimental 
feasibility. We indicate the free parameters with a star 
symbol. 


1. Nanodiamond 


* Radius nanosphere: i? = 15 nm. 

* Density of quantum emitters: pq = lA nm“^. This 
density is an upper value extracted from the exper¬ 
imental measurements in [48] , samples with smaller 
SiV densities would also fulfill Nr] ^ 1. 

• Number of quantum emitters: N = 1.98 x 10^. 

• Real part of the refractive index: Re[n] = 2.4 [77]. 
For the blackbody radiation, absorption and emis¬ 
sion, the leading term in the integral for moderate 
temperature {Tint, T^nv ^ 2000) arises from low en¬ 
ergy radiation (w < 2 x lO^^Hz). This constant 
value for the real part of the refractive index of 
diamond was consequently used. 


• Imaginary part of the refractive index in the spec¬ 
tral region 500-1100 nm: Im[n] = 10“®. This value 
is estimated from the values of extinction available 
in the literature (3.8 x 10“^ at 436 nm [77], < 
2x 10“® at 532 nm and < 9x 10“® at 1064 nm [78]). 

The spectral dependence of the imaginary part of 
the refraction used is given by: 


Im[n(a;)]=10 ®-|- 


+ j — + f2{x) + fs{x)\, 


where 


h{x) 

h{x) 

h{x) 


15 exp - 
1.5 exp 
0.35 exp 


2100- xV 
2 X 105 )\ ’ 
'3200-xV 
, 2 X 105 )\ ’ 
/4400-xY 
V 2 X 105 ) 


(Al) 


(A2) 

(A3) 

(A4) 


and X = w/(100 x 27rc). Such function provides 
a satisfactory fit of the experimentally measured 
values in the wavelength range 3-6.5^m [77]. In 
the absence of values for the extinction coefficient 
in the 6.5-20p,m, we assumed a value of 10“® which 
leads to a significant underestimation of the black 
body radiation of the nanodiamond. The complex 
dielectric function of diamond e (w) was obtained 
using this fit. 


• Transition frequency: a;o/27r = 4.01 x lO^'^ Hz. 
Its bulk temperature’s dependence is [50]: 


Ao(T) 

1 nm 


= 737-b 19.2 X 10 


-8 


T 

Ik 


2.78 


(A5) 


where Aq = 2t:c/ujo. 

Its bulk temperature’s dependence is [50]: 


r(r) 

277 X 109 Hz 


1 

9^ 


X _|_ 3.3e“55 nieV/{KhT) 


(A6) 


Its bulk temperature’s dependence is [50]: 

3 


liT) 


277 X 106 Hz 


= 16.39-b 1.9 X 10" 


I <"> 


3. Near-fleld trapping 

* Nano-diamond size: 30 nm. 

* Fiber diameter: 715 nm. 

* Detuning: A/27r = ± IxlO^^jjz. 

* Spatial decay of mode 1 (mode EH 21 ): 

A“^ =210 nm. 

-k Intensity for mode 1: Ii = 0.62mW//i^. 

* Spatial decay of mode 2 (mode HEn): 

A“9 =135 nm. 

* Intensity for mode 2: I 2 = 1.85mW//79. 

* Trap depth: 43 K. 

* Trapping distance from the fiber’s surface: 

Zt = 287 nm. 

* Nanodiamond internal temperature: 587 K. 


4. Microsphere cavity 

• Frequency of the cavity mode: uic = ujo + Ac where 
Ac = 1.4 X 10^5 112 

• Radius of the microsphere cavity: 25 /rm. 

• Cavity intrinsic radiative losses: Q~nd ~ xlO^®. 

• Cavity scattering losses: = 6.6 xlO^® [79]. 

• Cavity material losses: Q^at = 9 xlO^^ [79]. 

-k Cavity mode quality factor: Q = lO^^ < l/(Qj^s. + 

Qbulk T Qrad)' 
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• Decay rate of the cavity mode: 

K/27r = ujcl{‘2.TrQ) = 18.3 x 10^ Hz. 

• Cavity volume: 14 = 820 . 

• Decay of the evanescent field of the cavity mode: 
1/Ac = 283 nm. 

5. Resolved side-band regime 

• Distance between the nano-fiber and the cavity: 
D = 900 nm. 

X Detuning: A/27r = ± IxlO^^Hz. 

• Spatial decay of mode 1 (mode EH 21 ): 

=210 nm. 

• Intensity for mode 1: Ii = 0.62mW//r^. 

X Spatial decay of mode 2 (mode HEn): 

A“^ =135 nm. 

• Intensity for mode 2: I 2 = 1.85mW//r^. 

• Trap depth: 34 K. 

• Trapping distance from the cavity: 

Zt = 612 nm. 

• Nanodiamond internal temperature: 587 K. 

• Cavity loss rate due to the scattering from the emit¬ 
ters: Ksd K, = 0.34. 

• Mechanical decoherence rate due to photon scat¬ 
tering due to the emitters: r^/27r = 42.3 x 10^ Hz. 

• Mechanical decoherence rate due to photon scat¬ 
tering due to the nanosphere: r^/27r = 0.63 Hz. 

• Single-photon coupling vs. cavity loss rate: go/^ = 
1.2 X 10-4 

• Mechanical frequency: Wj/Itt = 1 x 10® Hz. 

• Mechanical frequency vs. cavity loss rate: wt/zc = 

4. 


6. Strong single-photon optomechanical 
cooperativity (FORT) 

* Distance between the nano-fiber and the cavity: 
D = 565 nm. 

* Far-of-resonance detuning: |A| = 1 x 10^® Hz. 


• Spatial decay of mode 2 (mode HEn): 

A“^ =81 nm. 

• Intensity for mode 2: I 2 = 1.43mW//r^. 

• Spatial decay of mode 3 (cavity mode): 

A-® =85 nm. 

• Intensity for mode 3: I 3 = 1.72mW///^. 

• Trap depth: 10.8 K. 

• Trap frequency: wt/(27r) = 32.7 x 10® Hz. 

• Trapping distance from the fiber’s surface: 

Zt = 270 nm. 

• Trapping distance from the microsphere cavity’s 
surface: z't = (565 — 270) nm. 

• Steady excited state population: pe = 3.2 x 10“®. 

• Nanodiamond internal temperature: T = 385 K. 

• Mechanical frequency vs. cavity loss rate: Wt/K = 

1 . 


• Single-photon coupling vs. cavity loss rate: gQln = 

2.8 X 10-4 

• Single-photon cooperativity: C = 1.2. 

• Cavity loss rate due to the scattering from the emit¬ 
ters: Ksc/k = 0.79. 

• Mechanical decoherence rate due to photon scat¬ 
tering from the emitters: r^/27r = 15.83 Hz. 

• Mechanical decoherence rate due to photon scatter¬ 
ing from the nanosphere: r^/27r = 76.5 x 10“® Hz. 

• Mechanical decoherence rate due to photon scat¬ 
tering from the black body radiation: r^^/27r = 

5.8 X 10"'^ Hz. 
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